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Abstract. This paper studies periodic traveling gravity waves at the free sur- 
face of water in a flow of constant vorticity over a flat bed. Using conformal 
mappings the free-boundary problem is transformed into a quasilinear pseudodif- 
ferential equation for a periodic function of one variable. The new formulation 
leads to a regularity result and, by use of bifurcation theory, to the existence of 
waves of small amplitude even in the presence of stagnation points in the flow. 
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1. Introduction 

The problem of spatially periodic traveling gravity water waves with constant vor- 
ticity 7 on a flow of finite depth can be formulated as the free-boundary problem of 
finding 

• a domain Q in (X, Y)-plane, whose boundary consists of the real axis 

(1.1a) B = {(I,0):l£l}, 

representing the flat impermeable water bed, and an a priori unknown curve 
expressed in parametric form as 

(Lib) S = {(u(s),v(s)) : s G R}, 

with 

(1.1c) u(s + L) = u(s) + L, v(s + L) = v(s) for all s G R, 

representing the free surface of the water, which is .L-periodic in the horizontal 
direction; 

• a function (X,Y) i— > t/j(X,Y) which is L-periodic in X throughout O, rep- 
resenting the stream function giving the velocity field (ipy, —ipx) in a frame 
moving at the constant wave speed, which satisfies the following equations 
and boundary conditions: 



(l.ld) Aip = -7 in SI, 

(Lie) ip = —m on B, 

(l.lf) ip = on S, 

(Llg) \V^\ 2 + 2gY = Q on S. 



Here g is the gravitational constant of acceleration, the constant m is the relative 
mass flux, while the constant Q is related to the hydraulic head (see the discussion 
in |12j). The level sets of ift are the streamlines, while a point where the gradient of 
ijj vanishes is called a stagnation point. 

We are interested in the existence, regularity and geometric properties of solutions 
of (|l.ip . The physical motivation for our study lies in the desire to understand 
wave-current interactions, since nonzero vorticity is the hallmark of a non- uniform 
underlying current \10\ 1 12 j . Even in the simplest case of a linear current, which 
corresponds to a constant vorticity 7 in (jl.ip . numerical studies [HI [22j [231 EZ] point 
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to the existence of large-amplitude waves with overhanging profiles and critical layers 
(regions in the fluid bounded by closed streamlines and containing stagnation points). 
This is in striking contrast to the extensively studied irrotational case [71 EH , for 
which any wave profile is necessarily the graph of a function and there can be no 
critical layers in the fluid |37[ I3U1 I3T] . The assumption of an underlying current with 
nonzero constant vorticity (or linear shear) provides us with the simplest case of a 
flow that is not irrotational and is attractive for analytical tractability. However, this 
setting is not a mere mathematical convenience, being physically relevant. Indeed, 
on areas of the continental shelf and in many coastal inlets the most significant 
currents are the tides and tidal flows are the most regular and predictable currents [19[ 
33j, considered to be two-dimensional flows of constant vorticity: negative constant 
vorticity is generated when the flow of the current is directed towards the shore 
(flood) and positive constant vorticity when the current is directed back out to sea 
(ebb) PS 05| . 

This paper addresses some of the mathematical challenges posed by the study of 
solutions of (|l.ip with overhanging profiles or critical layers. In all previous stud- 
ies of existence of solutions of (jl.ip . equivalent formulations of the problem over a 
fixed domain are essential. Bifurcation theory is then the most successful tool used 
for existence results. For irrotational flows (7 = 0, representing uniform flows with- 
out vorticity), a hodograph transformation maps the unknown domain occupied by 
the water into a fixed horizontal strip in a complex plane where the variable is the 
complex potential of the fluid flow. In these coordinates, the classical approach is 
to reformulate (jl.ip as a nonlinear singular integral equation for a function of one 
variable which gives the angle of inclination between the tangent to the free surface 
and the horizontal [U [2U1 A more recent approach El [31] (for the analogous 
problem of waves of infinite depth) uses a reformulation of the problem as a pseu- 
dodifferential equation for a function of one variable which gives the elevation of the 
free surface when the fluid domain is the conformal image of a strip [2J. 

In this paper we show that the second approach can be extended to the case of 
any constant vorticity 7. The basic idea is to seek the fluid domain as the conformal 
image of a strip, regardless of the unavailability of a complex potential. The new 
formulation is presented in Section 2. It imposes no restriction on the geometry 
of the free surface or on the streamline pattern, and in particular it can handle 
overhanging profiles and critical layers. In order to avoid distracting technicalities, 
we assume at the outset a modest amount of regularity of the free surface, namely 
C 1,a for some a 6 (0, 1), though arguing along the lines of |31[I41] one could verify the 
equivalence of the formulations also in the framework of weak solutions [43J . As the 
new formulation involves the periodic Dirichlet-Neumann operator and the periodic 
Hilbert transform associated to a strip, Section 3 studies some useful properties of 
these operators. Section 4 is devoted to regularity of solutions of the new equation, 
showing in particular that any C 1 '" free surface with no stagnation points is actually 
C°° . In Section 5 we prove the existence of waves of small amplitude using bifurcation 
theory. The analysis reveals the existence of such waves with a critical layer, a fact 
first proved rigorously only very recently in [33], though anticipated by a formal 
argument of Kelvin [34J, who presented the famous cat's eye picture of the streamline 
pattern. 

Before proceeding with the bulk of the paper we would like to comment upon the 
relation to some recent investigations on water waves with vorticity. For general vor- 
ticity but under the assumption that ip is monotone in the vertical direction, a hodo- 
graph transform (X, Y) 1— > (X, ip) was used in [12J to transform (jl.ip into a nonlinear 
boundary value problem for a quasilinear elliptic equation in a strip, which permits 
the investigation of waves of small and large amplitude (see also [91 [131 1121 S3]). Thus 
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the formulation in [T2] cannot describe neither critical layers nor overhanging pro- 
files. The recent paper [33] deals with waves of small amplitude in a flow of constant 
vorticity which could contain stagnation points. The approach of [33] cannot handle 
overhanging profiles. For a free surface of the form S = {(X,r](X)) : X G M}, it 
reformulates (jl.ip as a nonlinear boundary problem for an elliptic equation in a strip 



by performing the flattening change of variables (X, Y) i— >■ (X, ¥) = (X, 



The free surface is thus transformed into the horizontal line ¥ = 1 and, in the new 
variables, one has to solve a system of equations for a function of two variables and 
one of one variable. Our analysis of small amplitude waves is much simpler than 
that in [33]. The new formulation, which has an elegant structure, opens up the 
possibility of using global bifurcation theory to prove the existence of waves of large 
amplitude with critical layers and/or overhanging profiles, a prospect which seems 
daunting for the formulation in [33]. In addition, our approach has the advantage 
that in the reformulation we do not depart from harmonic function theory, and the 
fine structural properties associated with harmonic functions continue to be relevant. 
Such properties made it possible to describe the particle motion beneath irrotational 
waves [SI [T3] and the exploration of this aspect for flows with constant vorticity is 
now plausible. 



In this section we present the reformulation of the free-boundary problem (jl.ip 
as the quasilinear pseudodifferential equation (|2.9p for a periodic function of one 
variable. This involves the periodic Dirichlet-Neumann operator and the periodic 
Hilbert transform for a strip, so we start with a discussion of the definitions and 
basic properties of these operators. 

For any integer p > and a G (0, 1) we denote by C p,a the standard space of 
functions whose partial derivatives up to order p are Holder continuous with exponent 
a over their domain of definition. By C?_ we denote the set of functions of class C p ' a 
over any compact subset of their domain of definition. We will need the following 
local version of Privalov's Theorem [24\ Chapter V, §E, p. 100] (see \39\ Lemma 2.2] 
for the proof of a very similar result). A more refined, global version of Privalov's 
Theorem, related to but different from that in |24| . will be proved in Appendix B. 

Lemma 2.1. For any t > 0, we denote B+ = {(x,y) G M 2 : x 2 + y 2 <t 2 ,y> 0}. Let 
r > and Z + iW be a holomorphic function in where Z and W are real-valued 
functions. Suppose that W is continuous in B+ U {(x,0) : x G (— r, r)}, and let 



If w G C^'"((— r, r)) for some integer p > and a G (0, 1), then W and Z are of class 
C p,a in the closure of~Bf(xo) for every t G (0,r). 

For any d > 0, let TZd be the strip 



For any integer p > and a G (0, 1) we denote by C P £* the space of functions of 
one real variable which are 2-7r-periodic and of class C p ' a . For any w G C^, let 
W G C p ' a (TZii) be the unique solution of 




2. Reformulation of the free-boundary problem 



w (x) = W(x,Q) 



x G (— r, r). 



TZ d = {{x,y) el 2 : -d < y < 0}. 



(2.1) 



AW = 



in TZ d 



W(x,-d) = 0, 
W(x,0) = w(x) 
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The function (x,y) h-> W(x,y) is 2-7r-periodic in x throughout lZ d . For p > 1 integer, 
we define G d {w) by 

SdH(z) =W w (z,0), ieR. 
(Throughout the paper, the subscripts x, y, X, Y denote partial derivatives, and 
are the only subscripts with this meaning. ) Then g d (w) e Cf~ 1,Q! . The mapping 
w i — y G d (w) is called the periodic Dirichlet- Neumann operator for a strip, and is a 
(bounded) linear operator from to C\ v ,Q . Note in particular that if w is a 
constant function taking the value c, then 

W(x,y) = ~(y + d), (x,y)eK d , 

and hence 

(2.2) G d (c) = c/d. 

Suppose now that Z is a harmonic function in TZ d , uniquely determined up to a 
constant, such that Z + iW is holomorphic in lZ d . By Lemma 12.14 % G C p ' a (lZ d ). 
Note that the function 

(x, y) i ^ Z(x + 2vr, y) - y) 

is constant in 7^ as can be seen by taking its partial derivatives and using the 
Cauchy-Riemann equations. Let Kdlbe such that 

(2.3) Z{x + 2ir,y) -Z(x,y) = K, (x,y)eK d . 

For every i£l and y € (— d, 0), it follows using the Cauchy-Riemann equations that 

(2.4) —J W(t,y)dt = J W y (t,y)dt 

l'X + 2-K 

= / Z x (t,y)dt = Z(x + 27r,y)- Z(x,y) = K. 

J X 

Integrating the above relation from y = — d to y = 0, and taking into account the 
boundary conditions satisfied by W in (|2.ip . it follows that 

(2.5) K- 27rM 



d 1 

where [w] denotes the average of w over one period. Hence K = if and only if 
[w] =0. For any integer p > and a € (0, 1), let us denote by Cf^ a the class of 
functions in C^" which have zero mean over one period. 

Let w G C°2^ . Then the function (x,y) h-> Z(x,y) is 27r-periodic in x throughout 
lZ d . We normalize the constant in the definition of Z by the requirement that x i— > 
Z(x,0), which is a 2-7r-periodic function, has zero mean over one period. We define 
C d {w) by 

C d (w){x) = Z(x,0), i£R. 

The mapping w t— >■ C d (w) is called the periodic Hilbert transform for a strip, and is a 
(bounded) linear operator from into itself. It is easy to see from the definitions 
and the Cauchy-Riemann equations that, when w E Cf^o f° r P — 1 integer, 

&(«;) = (CdCw))' = QK). 

A more detailed account of the properties of the operators Q d and C d is given in 
Section 3. Note that the space of 2-7r-periodic harmonic/holomorphic functions in 
a strip can be identified with the space of harmonic/holomorphic functions in an 
annulus, by means of an explicit conformal mapping, and thus the operators Q d and 
C d can be equivalently defined in the setting of an annulus. 
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Suppose now that w £ C^*. By (|2.2p . (j2.3p and (|2.5p . one can write 
Z + iW = ^(x + i{y + d)) + Z + iW mTZ d , 

where Zo,Wq £ C p ' a (1Zd) are 27r-periodic in x throughout 1Z d and 

(2.6) AW = in K d , 

W (x,-d) = 0, x£R, 

W (x,0) = w(x) - [w], xel. 

We normalize the constant in the definition of Z by the requirement that x t— >■ Zq(x, 0), 
which is a 2-7r-periodic function, has zero mean over one period. It follows that 

(2.7) Z(x,0) = ^-x + C d (w - H). 

It also follows, using the Cauchy-Riemann equations, that, when p > 1, 

(2.8) g d ( w ) = M + - H)y = M + c d K). 

Throughout the paper we are interested in solutions (fi, ?/>) of the water-wave prob- 
lem (jl.ip of class C 1,a , for some a £ (0, 1), by which which we mean that S has a 
parametrization (jl.lbp with u, v functions of class C 1,a , such that (jl.lcp holds and 

u'{sf + v'{sf ^ for all s £ R, 

while ^ £ C°°(r2)nC 1 ' a (0). The main result of this section is that the free-boundary 
problem (jl.ip is then equivalent to the problem of finding a positive number h and a 
function v £ C 2 £* which satisfy the following: 

(2.9a) + j(g kh (v 2 /2) - vQ kh {v)) f = (Q - 2gv) [v' 2 + Q kh {vf) , 

(2.9b) [v] = h, 

(2.9c) v{x) > for all x £ M, 

(2.9d) the mapping x i-)- ^— +Ckh{v — h)(x),v(x) S j is injective on R, 

(2.9e) v'{xf + Qkh{v){xf ^ for all x £ R, 

where 

(2.10) L = 2Tr/k, k>0. 

(Thus 2-/r-periodicity in x becomes L-periodicity in s = x/k.) More precisely, we have 
the following. 

Theorem 2.2. Let be a solution of (jl.ip of class C l,a . Then there exist a 

positive number h, a function v £ C 2 ^ and a constant a £ R such that (|2.9p holds 
and 

(2.11) s = {(a+^ + C kh (v - h)(x),v(xfj : x £ R} . 

Conversely, let h > and v £ be such that (|2.9p holds, and let a £ R 6e arbitrary. 
Let S be defined by (|2.1ip . and Q 6e t/ie domain whose boundary consists of S and 
the real axis B. Then there exists a function tp in fi such that (O,^) is a solution of 
(JED o/c/ass C7 1 '". 

The presence of the constant a in (|2.1ip is due to the invar iance of problem (|1.1|) 
to horizontal translations. We now explain the meaning of the positive constant h 
which arises in (|2.9|) . A domain 17 contained in the upper half of the (X, Y)-plane is 
called an L-periodic strip-like domain if its boundary consists of the real axis B and 
a curve S described in parametric form by (jl.lbp such that (jl.lcp holds. For any 
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such domain, we define its conformal mean depth as the unique positive number h 
such that there exists a conformal mapping U + iV from the strip IZh onto fi, which 
admits an extension as a homeomorphism between the closures of these domains, 
with {(x, 0) : x G E} being mapped onto S and {(x, —h) : x G R} being mapped onto 
,6, and such that 

(2.12) U(x + L,y) = U(x,y) + L, V(x + L, y) = V(x, y), (x,y)eTZ h . 

The existence of a unique number h with this property is proved in Appendix A. (The 
notion of conformal mean depth is reminiscent of that of the conformal modulus for 
doubly connected domains, defined as the unique number R G (0, 1) such that the 
domain is conformally equivalent to the annulus {(x,y) G R 2 : R 2 < x 2 + y 2 < 1}, 
see \26\ I30|.) Note that, for any horizontal strip, its conformal mean depth coincides 
with its usual depth (height). Also, as we show in Appendix A, the conformal mean 
depth coincides with the 'mean depth' used in \\\ in the study of irrotational flows. 
The proof of Theorem 12.21 reveals that h in (|2.9f) is exactly the conformal mean depth 
of O. 

Proof of Theorem \2.2l Let be a solution of (jl.ip of class C 1,a . Let h be the 

conformal mean depth of Q, and let U + iV be the associated conformal mapping. 
As we show in Appendix A, it is a consequence of the Kellogg- Warschawski Theorem 
that U,V€ C l ^ a (jZh) and 

U 2 (x, 0) + V 2 (x, 0) / for all x G M. 

Consider the mapping U + iV : IZkh ^ given by 

(2.13) U(x,y) = U(x/k,y/k), V(x, y) = V(x/k, y/k), (x,y)eTl kh , 

where k is given by (|2.10p . Then U + iV is a conformal mapping from IZkh onto 
17, which admits an extension as a homeomorphism between the closures of these 
domains, with {(x,0) : x G M} being mapped onto S and {(x,—kh) : x G M} being 
mapped onto i3, and such that 

2ir 

(2.14) [/(x + 27r,y) = [/(x,y) + — , V{x + 2vr, y) = V(x, y), (x,y) eK kh . 
Moreover, U, V G C x ' a (^) and 

(2.15) C/J(x, 0) + V 2 (x, 0) / for all x G R, 
while the conformal mapping properties of U + iV imply that 

(2.16) V(x,-h) = Q, 

(2.17) the mapping x h-> (U(x,0),V(x,0)) is injective on R, 

(2.18) S = {(£7(x,0),T/(x,0)) : x G R}. 
Let 

(2.19) u(x) = T/(x, 0) for all x G R. 

Then v G . Since V is harmonic in IZkh and satisfies (|2,16p and (|2.19p , we deduce, 
by comparing (pJ4l with ([23]), §Mk and ([23]), that 

(2.20) [v] = h. 

It follows, by using (|2.7p . that ()2.17p and ()2.18p can be rewritten as 

(2.21) the mapping x t-t ^— + Cfc^(w — h)(x), v(x) S j is injective on R, 

(2.22) 5 = { (a + | + C fe7l (<; - h)(x),v(x^J : x G r} , 
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for some aeK. Also, using the Cauchy-Riemann equations one can rewrite (|2.15|) as 

(2.23) v'(x) 2 + Q kh {v){x) 2 / for all x G R. 

Observe also that, since S is contained in the upper half-plane, it follows that 

(2.24) v{x) > for all i£l, 

Suppose now that ip satisfies (jl.ip . and let £ : TZ k h — > R be given by 

(2.25) £(x, y) = HU(x, y), V(x, y)), (x, y) G TZ kh . 
Note that (jl.ldp can be rewritten as 

(X, Y) i — y ip(X, Y) + t^Y 2 is a harmonic function in f2. 

Hence, since harmonic functions are invariant under conformal mappings, it follows 
that 

(2.26a) £ + is harmonic in TZ k h- 

It is immediate from (jl.lep and (|l.lf|) that 

(2.26b) £(x, -kh) = -m, x G R, 

(2.26c) £(je,0) = 0, 

By using the chain rule and the Cauchy-Riemann equations, we easily obtain that 

£ + f y = Wx(U, V) + j&(U, V))(V 2 + V y 2 ) in TZTh- 
It follows from Ql.lgP , taking into account (|2.15p . that 
(2.26d) il + il = (Q - 2gV)(V 2 + V 2 ) at (x, 0) for all i£i 

Let C : ^-fch -> K be given by 

(2.27) C = z + m + lv 2 . 

Then ()2.26p can be equivalently rewritten as 

(2.28a) AC = in TZ kh , 

(2.28b) ((x, -kh) = for all x G M, 

(2.28c) C(x, 0) = m + |w 2 (x) for all x G R, 

(2.28d) (£, - ~/VV y ) 2 = {Q- 2gV)(V 2 + V 2 ) at (x, 0) for all x G R. 

Equation (|2.28|) can be conveniently expressed by means of the Dirichlet-Neumann 
operator as 

(2.29) {^+l{Gkh{v 2 /2)-vg kh (v)^y = (Q- 2gv) (y 12 + Q kh {v) 2 ^ , 

where we have used the linearity of and the fact that its action on constant 
functions is given by l|Z3]>.Bv gathering flZZSD , dZZID , (|2I55|> . ([ZMD and jZZSH , we 
obtain ([23]) . 

Conversely, suppose that the positive number /i and the function v G C 2 ^ satisfy 
(I2.9p . We now show how one can construct a solution of (jl.ip by reversing the process 
which led from (jl.ip to (j2.9|) . Let V be the harmonic function in TZkh which satisfies 
(12T6D and (127T9L and let U : TZ kh -> R be such that £7 + iV is holomorphic. By 
LemmaEH U, V € C 1 '"(^^). Condition (jHEj) ensures that (jgHD holds. Condition 
(|2~9d|) shows that the curve 5 given by (|2.18p . which can also be expressed as (|2.22|) . 
is non-self-intersecting, and (|2.9c[) shows that it is contained in the upper half-plane. 
(The fact that U is determined by V only up to an additive constant leads to the 
presence of an arbitrary constant a G R in (|2.22p .) Let Q be the domain whose 
boundary consists of 5 and B. An application of the Darboux-Picard Theorem [U 
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Corollary 9.16, p. 310], see [41 [ Proof of Theorem 3.4] for details, shows that U+iV is a 
conformal mapping from TZ^h onto $7, which admits an extension as a homeomorphism 
between the closures of these domains, with {(cc,0) : x S M} being mapped onto 
S and {(x,—kh) : x G M} being mapped onto B. It follows from (|2.14|) that $7 
is a L-periodic strip-like domain, where L is given by (|2.10p . The domain $7 has 
conformal mean depth h, since the mapping U + iV : TZh — > J7, where £7, V are given 
by (|2.13p . has all the required properties. Note also that (|2.9ep shows that S is a 
C 1,a curve. Let us define £ as the unique solution of (|2.28a|) - (|2.28cp . Obviously, 
C G C 1 '°(^) n C°°(TZ kh ). Let us define £ by (|227)) . and then V by (T2T25D . We 
obtain that V G C 1 '"^ n C°°(Q) satisfies (fTl^ - (fl~Tl . Finally, since (f2T9aD holds, 
we obtain that ip satisfies Ql.lgP - This completes the proof. □ 



free surface X = u (s) , Y = v (s) 



y = 

¥(X,Y) ..... U + iV 



% (x,y) 



flat bed Y = y = - kh 



Figure 1 . The conformal parametrization of the fluid domain . 

For the sake of completeness, notice that differentiating (|2.25p . solving the resulting 
linear system for ipx and ipy and taking subsequently (|2.27j) into accoount, we find 

that the velocity field (ipy, — ipx) at the location (X,Y) = (lJ(x,y), V(x,y) S j £ $7, 

where (x,y) G Tlkh-, can be expressed as 

VxCx ~t~ VyCy \t VxCy VyCx 

V 2 + V 2 7 V 2 + V 2 

r x 1 " y v x 1 * y 

in terms of C,(x, y) and of the conformal map U + iV from IZkh to fL 

3. On periodic harmonic functions in a strip 

This section contains a more detailed investigation of the operators and 
for d > 0, in particular their representation as singular integrals and some of its 
consequences. These results are important in the study, carried out in Sections 4 and 
5, of the regularity and local bifurcation of solutions of (|2.9p . 

Let be the space of 2-7r-periodic locally square integrable functions of one real 

1 2 

variable, and W 2 ^ be the space of 2-7r-periodic locally absolutely continuous functions 
whose weak derivatives are in L 2 n . By and W 2 ^ a be denote the subspaces of L 2 W 

and whose elements have zero mean over one period. Every function w G L 2 n 
has a Fourier series expansion 

oo oo 

(3.1) w = [uu] + a n cos(rax) + b n sin(nx), 

n=l n=l 
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and, if we denote 



\ 1/2 

12 



w \\ = [ W l 

V n=l / 

2 



then, by Parseval's Identity, | ■ || is an equivalent norm to the standard one in L^r- 
For any d > and w G L 2 .^, the function W : TZd — > R given by 



a sinh ma ^-^ sinh na 

n=l n=l 

is the unique solution of (|2.ip with the third condition there being satisfied not 
necessarily in the classical sense, but rather as 

(3.2) lim ||W(-,y) - toll = 0. 

When w G W^, we define Qd{w) as the unique function in such that 
(3.3) 

namely 

r 1 OO OO 

(3.4) Qd{w) = —j- + nCLn c °th(n(f) cos(nx) + nb n coth(nd) sin(nx). 

n=l n=l 

The mapping is the periodic Dirichlet-Neumann operator for a strip, a bounded 
linear operator from W 2 ^ into L\^. 

Suppose that w G Q . Then any harmonic function Z in TZd such that Z + iW 
is holomorphic is given by 

Z(x, y) = C + V] Q ^ c °sh(n(y + dj sin ( nx ) _ cosh n(y + rf) cos ( nx ^ 

smh(nd) smhnd 

n=l n=l 

where C is a constant. We choose C = and define Cd{w) as the unique function in 
o such that 

(3.5) lim||Z(.,j/)-C d (ti;)||=0, 
namely 

OO OO 

(3.6) = a n coth(nd) sin(nx) — 6 n coth(nd) cos(nx). 

n=l n=l 

The mapping is the periodic Hilbert transform for a strip, a bounded linear oper- 
ator from into itself. Moreover, one can easily see that Cd is a bijection from 
L^tt o onto itself. Let us denote by C^ 1 its inverse. Observe that, by formally setting 
d = oo in (|3.6p . we obtain the familiar periodic Hilbert transform [71 131} 38j. 

OO OO 

(3.7) C{w) = a n sin(nx) — b n cos(nx). 

n=l n=l 

for all w G with the Fourier expansion (|3.ip . The operator C is a bijection from 
L^tt o onto itself, and C _1 = — C. 

The operator C has a pointwise almost everywhere representation as a singular 
integral 

(3.8) C(w)(t) = ^PV J* cot (tli\ w (s)ds, 

where PV denotes a principal value integral [32]. The representation (j3.8|) is useful 
in proving Privalov's Theorem [23]: for every a G (0, 1) the operator C is a bounded 



(3.3) lim\\W y (;y)-g d (w)\\=0, 

y/*0 
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linear operator from C^f 1 into itself. Moreover, (|3.8|) is also instrumental in proving 
(see Lemma 10.5.3, p. 148]) that the operator 

(3.9) w i-> Q(w) = w C(w') - C(ww') 

maps functions w E C^f 1 into functions in for any 5 E (0, a). We now prove that 
these results, and the corresponding ones for functions of higher regularity, extend to 
the case of Cd, for any d > 0, by exhibiting an analogous representation of Cd as a 
singular integral. Let us denote by Qd the mapping 

(3.10) w^Q d (w)=wC d (w')-C d (ww'). 

Lemma 3.1. For any d > 0, p > integer and a E (0, 1), Cd is a bounded linear 
operator from into itself. Moreover, C^ 1 is also a bounded linear operator from 

into itself. 



or 



Lemma 3.2. If w E Cf^* with p > 1 integer and a E (0, 1), then Qd(w) E C^f /< 
any 5 E (0, a). 

Observe that the operators C and C d can be extended by (complex-)linearity to 
complex- valued functions in L\ v , and are characterized by their action on the 
trigonometric system {e mt } nG z\{o} : 

(3.11) C(e int ) = -i sgn(n) e int , n E Z \ {0}, 
and 

(3.12) C d (e int ) = -i coth(nd) e int , n£Z\ {0}. 
Let us write 

(3.13) C d = C + K d . 

Then the operator ICd corresponds to the Fourier multiplier operator on L\ v Q : 

(3.14) E c ne mt }^{ £ -isgn(n)A„ Cri e m '} 

n€Z\{0} n£Z\{0} 



with 

X — 

g 2|n|d _ J 



A„ = -^t-t-; r, |n| > 1. 



Notice that for any n > 1, 

oo 

^m! " (p+1)! ^" 



yields 



so that 



n p ^ i (p + 1; 



^ e 2nd _ 1 - n (2d)P+l ' P - 
n€Z\{0} v ; neZ\{0} 



Let Atrf E Lf^.o be given by 

oo 

(3.15) K d (t)= Yl -isgn(n) X n e int = ^2A n sin(nt), tel. 

n£Z\{0} n=l 
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Since the Fourier coefficients of n d decays faster than any power of \n\, n d is a function 
of class C°° (see [17]). It follows from (|3.14p that JC d (w) is the periodic convolution 
of w with the smooth function K d given by (|3.15p . namely 

i r 

(3.16) JC d (w)(t) = — / K d (t - s)w(s)ds, tel. 

Proof of Lemma V3.1\ It is easy to see that JC d is a bounded linear operator from O^ 
into itself. The operator C is also a bounded linear operator from into itself, 

by Privalov's Theorem and the fact that C commutes with differentiation on C 2 £* - 
The required result for C d follows from (|3.13p . The corresponding result for C^ 1 is 
obtained by the same argument, after observing that one can write 

(3.17) Cj 1 = -C + K d , 

where fC d is given by periodic convolution with a smooth function. 

□ 

Proof of Lemma \3.2l Using (|3,10p and (|3.16p , we write 

(3.18) Q d (w) = Q(w) + {w K d {w') - K d {ww')}. 

Since n d is smooth, it is easy to see from (|3.16p that JC d (w') and K, d (ww') are both 
smooth functions. The required result follows then from (|3.18p . provided that one 
can show that Q maps C P ^ into for any 5 G (0, a). This fact has been proved in 
[7J for p = 1, and in [28] for p > 2, the proof in [28] being in fact in a more general 
setting. 

For the sake of completeness, we now give a self-contained proof of the fact that, 
for any p > 1 integer, Q maps into for any 5 G (0,a). Our proof relies 
on a commutator estimate proved in Appendix B, which is more general and has a 
more transparent proof than the results in [28]. Let w G C^f , with p > 1 integer and 
a G (0,1). Since G , we may interchange the operators C and differentiation 
to infer from Leibniz's rule that 

(wC(w {p) ) -C(ww {p) )^j - SfP^fwCiw') - C(W)) G Cl; . 
Thus it suffices to show that for any 5 G (0, a) we have 

(3.19) (wC(w^) - C(ura/ p ))) G G^. 

However, estimate (|3.19p is an immediate consequence of Lemma B.l in Appendix B. 
This completes the proof. □ 

Remark 3.3. Since for a G (0, 1) the function 

W : [-7T, 7T] -> M, W(X) = - 7T Q 

extends by periodicity to a function to G with wC(w') — C{ww') G" C 2 ^* (see 
the detailed estimates in [28J), in view of (|3.16|) we see that the result of Lemma 3.2 
cannot be improved to accommodate 5 = a. 

4. Regularity 

The main result of this section concerns regularity of solutions of (|2.9p . 
Theorem 4.1. Let h > and v G be a solution of (j2.9j) such that 
(4.1) Q-2 5 u(x)>0 fordlx£R. 

Then veC£. 
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Remark 4.2. Theorem 14.11 implies that, is is a solution of (jl.ip of class C ,a 

with no stagnation points on the free surface, then S is a C°° curve and ip G C°°(f2). 
Indeed, by Theorem 12 .21 any such solution gives rise to a solution of (|2.9p . the absence 
of stagnation points on S being equivalent to (|4.ip . By Theorem 14. 1\ v G and, 
by Lemma IBTTl Ckh{v — h) £ C^. Hence S, being given by (|2.1ip . is a C°° curve, and 
the fact that tp G C°°(f2) follows from standard elliptic regularity theory 



Remark 4.3. Our proof of Theorem 14.11 is based on harmonic analysis estimates 
and is relatively simple. Using different methods which are somewhat more involved, 
the stronger conclusion that s i— > v{s) is real-analytic on R can be obtained. Indeed, 
once the C 2,a regularity of S and of ip in O is known (for example from Theorem 
4.1), the approach of Kinderlehrer, Nirenberg and Spruck [21j is applicable, as first 
observed by Constantin and Escher [11], showing that the curve S is real-analytic. 
Then a standard result of regularity of conformal mappings up to the boundary [29] 
shows that U + iV has a holomorphic extension to a neighbourhood of the real axis, 
implying in particular that s t— >■ v (s) is real-analytic on KL An alternative proof of 
this fact, based on ordinary differential equations in the complex domain in the spirit 
of Lewy's theorem [25], has been pointed out to us by the anonymous referee of this 
paper. 



Proof of Theorem \4-l\ Let U, V be as in the proof of the second part of Theorem I2.2| 
with U, V G C^ a {U kh ). As noted there, U + iV is a conformal mapping from IZkh 
onto fi, and therefore 

U x + iV x 7^ in TZ kh . 

Since V > in TZkh an d V(x, —kh) = for all x E K, it follows from Hopf boundary- 
point lemma [18] that 

V y (x, -kh) > for all x£l, 
Taking also into account ()2.9ep . we have thus obtained that 



(4.2) U x + iV x ^0 mK kh . 
One can therefore write 

where 6 G C°' a (7l^) is such that 

log {Ul + V 2 ) 1/2 + iO is holomorphic in TZ kh , 
9{x, -kh) = for all x G E. 

Let 

(4.3) 6 (x) = 6(x, 0), u(x) = U(x, 0), 
Then 

(4.4) C kh (6 ) = log (u' 2 + t/ 2 ) 1 /* - [log (u' 2 + v' 2 ) 1 ' 2 ), 
so that 

(4.5) 9 = Cll (log (u' 2 + ^' 2 ) 1 / 2 - [log (u' 2 + v' 2 



Since ()4.ip holds, equation ()2.9ap can be written as 
(4.6) {u^ + v' 2 ) 1 ' 2 



^ + l{Gkh(v 2 /2)-vg kh (v)} 



(Q - 2gv)V 2 

Notice that in view of (J23D, (f2T9b]l and (l3TT0]l . we have 



(4.7) S fe/ > 2 /2) - w g fcfc («) = |J + C kh (vv') - j -vC kh (v') = t±-l- Q kh ( v ). 

Ikh k 2kn k 
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Since v G C^, Lemma [3721 ensures that Qkh(v) G and from (|4.6p and (|4.7p 

we infer that (u' 2 + v' 2 ) G Ci^ 2 . Using this in (|4.5p . Lemma 13.11 ensures now that 
6 G Cl' n a/2 . But 

|y = (««+ i^)V2 S i n O) 
so that u,i)£ • Implementing this procedure p times we obtain that 

u,t,€Ci^^, P >1. 
Consequently u G C^. This completes the proof. □ 

5. Local bifurcation 

In this section we prove the existence of solutions of (|2.9p . Because of (|2.9bp , it is 
natural to put 

(5.1) v = w + h, 

and rewrite (|2.9p . taking into account (j2.8|) . as 

f m ( \w 2 } w h „ . „ . 

(5.2a) \kh +7 \2kh ~I~2k + C - (W) " wCkh{w) 

= {Q- 2gh - 2gw) ju/ 2 + + Ckh(w' 

(5.2b) M = 0, 

(5.2c) it; (a;) > -h for all i£K, 

(5. 2d) the mapping x h-> ^— + Ckh{w){x), w(x) + /i^ is injective on R, 

(5.2e) w/(x) 2 + Q + Cfefc(«/)(a:)^ ^ for all x G R, 

We prove, for each 7 G R, > and h > fixed, and regarding m and Q as 
parameters, the existence of solutions w G C 2 ' Q of (|5.2|) . Note that, if 

(5.3) - + Cm(OOc) > for all xeR, 

k 

then (]5.2dp and ()5.2ep are automatically satisfied and, moreover, the corresponding 
free boundary S is the graph of a function. Motivated by the remark that w suf- 
ficiently small in ensures the validity of ()5.2cj) and (|5.3p . we concentrate on 
proving the existence of small-amplitude solutions w G C 2 !^ of (|5.2ap . Of course, 
any study of large-amplitude solutions of (|5.2ap - (j5.2bp will need to discuss also the 
validity of ()5.2cp - (|5.2ep . which is necessary for the construction of solutions of (jl.ip . 
Note that w = G C 2 £* is a solution of (|5.2a|) if and only if 

(5.4) Q ^ h+ (^-±)\ 
This suggests setting 

m 7/1 

(5.5a) ~h~\' 

(5.5b) , 1 = Q-2gh-('^-^)\ 
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The mapping (m, Q) i— > (A, //) is a bijection from R 2 onto itself. Equations (|5.2aj) - 
()5.2bj) can be rewritten as 

(5.6) |^ + 7 ~j + C kh (ww') - wC kh (w') 

= (A 2 + n - 2gw) iw' 2 + Q + C k h(w')j j , 

with w G C^q, fJL G R and A € R. Note that w = 6 C^tTo an d A* = is a solution 
of flSJl]) for every A G R. 

To prove the existence of solutions of (|5.6p we apply the Crandall-Rabinowitz 
theorem [15] on bifurcation from a simple eigenvalue. For a linear operator £ between 
two Banach spaces, let us denote by M(C) its null space and by TZ(C) its range. 

Theorem 5.1 (Local bifurcation theorem [15j). Let X and Y be Banach spaces, 
I an open interval in R containing A*, and F£C(JxX,Y), Suppose that 

(i) F(A,0) = /or a// A G I; 

(ii) d\F , d u F , and d 2 u F exist and are continuous; 

(Hi) M(d u F(X* , 0)) and Y /lZ(d u F(X* , 0)) are one- dimensional, with the nullspace 
generated by u* ; 

(iv) the transversality condition d 2 U F(X*, 0) (1, u*) G" 1Z(d u F(X* , 0)) holds. 
Then there exists a continuous local bifurcation curve {(A(s),w(s)) : |s| < e} with 
e > sufficiently small such that (A(0),u(0)) = (A*,0) and 

{(X,u) GO: u ^ 0, F(X,u) = 0} = {(A(s),u(s)) : < |s| < s} 

/or some neighborhood O of (A*,0) GJxX. Moreover, we have 

u(s) = su* + o(s) m X, \s\ < e, 

and if d 2 F is also continuous, then the curve is of class C 1 , while for F of class C k 
(k > 2) or real- analytic, s i— > u(s) is of class C k ~ 1 , respectively real- analytic. 

Let us denote, for every integer p > 0, 

(5.7) C*£ = {/ G Cg? : f(x) = f(-x) for all x G R}, 

(5.8) e^ 0>e = {/ G e 2 ^ : /(*) = /(-*) for all x G R}. 
To apply the local bifurcation theorem to (|5.6p . let 

(5.9) x = r x eg*- y = c*£, 

for some p > 0. Then we can write (|5.6p as -F(A, (fi,w)) = with FilxX-^Y 
real-analytic given by 

(5.10) F(\,(ji,w)) = 7 2 fc fe?l (W)-wC fc / l K)-| + ^) 

2A7 / . , n I l\ W 



+ 



; (e kh { ww ') - w c kh ( w ') -j + 2lh 



+ (2gw - n) + w' 2 + ^ C kh (w') + ( C/,./, ( 

-A 2 ^' 2 + ^C fc/t K) + (c fc/l K)) : 

As noted earlier, F(A, (0,0)) = for all A G R. Taking into account Lemma I3TT1 we 
easily compute 

(5.11) d M F(X, (0, 0)) (/, u) = ^ (( 5 - A 7 )/ - X 2 kC kh (f')) - ^, (f, v) G X. 
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It follows from the representation (|3,6|) and Lemma 13.11 that the bounded linear 
operator dr^w) F(X, (0,0)) : X — > Y is invertible whenever A does not satisfy 

(5.12) A 2 n/ccoth(n£;/i) = g - A7, 

for any integer n > 1. Hence all potential bifurcation points for (|5.6p are to be found 
among the solutions of (j5. 12[) for some integer n > 1. 

Suppose now that A* is a solution of (|5.12p for some integer n > 1. (Observe 
that (|5.12p has exactly two different solutions for each n > 1.) Using again the 
representation (|3.6p and Lemma [3.1 1 it follows that M(d^^F(X* , (0, 0))) is one- 
dimensional and generated by (0, w*) £ X, where w*(x) = cos(nx) for all x € R, 
while TZ(d^ w }F(\* , (0, 0))) is the closed subspace of Y formed by the functions / € Y 
satisfying 

t"IT 

/(x) cos(nx) dx = 0, 



so that Y /TZ(d^^ w )F(X* , (0,0))) is the one-dimensional subspace of Y generated by 
the function w*(x) = cos(nx). Using (|5.1ip . we now compute 

d t M H>f, (0,0)) (1, (0,tO) = ^(-7-2A*nfccoth(nA:/ l )) w* £* TZ(d M F(X*, (0,0))) 
since, using (|5.12|) . we have 

—7 — 2X*nk coih(nkh) = — A* [nk coth(nkh) + ) ^ 0- 

Therefore, by the local bifurcation theorem, the solutions of (|5.12p for any integer 
n > 1 are bifurcation points. Motivated by the quest for solutions of (|5.6p of minimal 
period 2ir, we take n = 1 in ()5.12p . obtaining the bifurcation values 



7tanh(/c/i) / 7 2 tanh 2 (/c/i) tanh(A;/i) 

(5 - 13) A± = 2^ ± V 4P + 9 ^— 

The corresponding values of m are obtained from (|5.5a|) . 



7/1 2 7/itanh(A;/i) / 7 2 tanh 2 (A;/i) tanh(/c/i) 

(5 ' 14) m ± = — ~ 2k ~ ±h f J 

The existence of water waves of small amplitude is now immediate. 

Theorem 5.2. Given h > 0, k > 0, 7GR anc? m £ 1R i/iere exist laminar flow£\ with 
a flat free surface in water of depth h, of constant vorticity 7 and relative mass flux m. 
Moreover, the values m± of the flux given by (|5.14j) trigger the appearance of periodic 
steady waves of small amplitude, with period 2ir/k and conformal mean depth h, which 
have a smooth profile with one crest and one trough per period, monotone between 
consecutive crests and troughs and symmetric about any crest line. The laminar flows 
of flux m± are exactly those with horizontal speeds at the flat free surface X± given 
by (15TT3L 

Proof of Theorem \5.2[ The function w = satisfies (|5.2p for any value of m £ R., 
provided that Q is given by (|5.4p . These solutions correspond to laminar flows in the 
fluid domain bounded below by the rigid bed B and above by the free surface Y = h, 
with stream function 

^(X,Y) = -lY 2 +(B + ^\ Y -m, X € R,0 < Y < h, 
and velocity field 

(5.15) (th,-tpx)= (-^ y + X + ^'°)' XeR,Q<Y<h. 



4n the sense that the water flows in layers parallel to the flat bed (with zero vertical fluid velocity). 



16 



ADRIAN CONSTANTIN AND EUGEN VARVARUCA 



Observe that for these flows the horizontal velocity at the free surface is indeed A 
given by (|5.5ap . 

Consider now the local bifurcation curve 

{(A(s), (0 + o(s), s cosO) + o(a))) : |s| < e} C R x X 

of solutions of (|5.6p . issuing from the points (A±, (0,0)), with A± given by (|5.13p . If 
e > is chosen small enough, then conditions (|5.2c|) and (|5.3p are satisfied, where we 
have used Lemma 13. 11 Therefore the corresponding non-flat free surface S given by 
(gT[]) with 

v = w + h, 

is the graph of a smooth function, symmetric with respect to the points corresponding 
to x = mr, n G Z. If we choose p in (|5.9p such that p > 1, then, since 

(5.16) w(x; s) = s cos(x) + o(s) in ' a , 
we can ensure that 

sw'(x; s) < for all x G (0, 7r), < |s| < e, 

by choosing e > suitably small. Together with the evenness of x i— > w(x;s), this 
proves the S has one crest and one trough per minimal period and is monotone 
between consecutive crests and troughs. 

□ 

Let us consider the family of laminar flows from which bifurcation of small-amplitude 
periodic waves was shown to occur. It follows from (|5.15p that they have the form 

(5.17) W>y,-Vx) = U± + l(h-Y),0\, X G M, < Y < h. 

where A± is given by ()5.13p . The formula (|5.13p . giving the speed A± at the free 
surface in terms of the depth h, period 2n/k and vorticity 7, is called the dispersion 
relation. As far as we are aware, its first occurrences in the literature are in [35] and 
where it was obtained by formal arguments. 
A remarkable feature of some of these flows is that they contain stagnation points. 
Since the flows are laminar, stagnation points, if present, form horizontal lines. An 
elementary analysis, similar to that in [35], shows that small-amplitude waves bifur- 
cating from these flows have a critical layer. The streamline pattern is considerably 
different in the case when stagnation points are present compared to the case when 
they are not. A typical case of a flow without stagnation points is depicted in Figure 
2 (see the discussion in [12]). 



free surface 




flat bed 



Figure 2 . Flow pattern (streamlines) without stagnation points 
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A picture of the streamline pattern if the flow presents stagnation points is provided 
in Figure 3: the 'cat's eye' flow pattern of Thomson (Lord Kelvin) [M] is highly 
typical for this situation (see the discussion in [S]). 



free surface 




flat bed 

Figure 3 . Flow pattern (streamlines) with stagnation points : 
Kelvin's 'cat's eye' flow pattern 



Notice that A+ > and A_ < irrespective of the values of h, k and 7, so stagnation 
points can never occur at the free surface of the bifurcation-inducing laminar flows. 
This implies that in the irrotational case (7 = 0) stagnation points do not occur 
anywhere in these flows. In the general case stagnation points occur if and only if 
the equation 

A± + 7 (/i-y) = 

has a solution Y in [0, h]. It is easy to see that this is the case if and only if 

(5.18) A±(A± + 7/1) < 0. 

Suppose first that 7 > 0. Then A + + 7/1 > 0, so the flow corresponding to A + never 
contains stagnation points. The flow corresponding to A_ contains stagnation points 
if and only if A_ + jh > 0, which is equivalent to 

(5.19) t&nHkh) < 72/1 



kh g + 7 2 /i 

For fixed h, the left-hand side of (|5.19p is monotone as a function of k, with limit 1 as 
k — > and limit as k — > 00, while the right-hand side is monotone as a function of 
7, with limit as 7 — > and limit 1 as 7 — > 00. Hence, for fixed h and 7, there exists 
a unique k* > for which equality holds in (|5.19p . which corresponds to a laminar 
flow stagnant at the bottom, and for k < k* the flow does not have stagnation points, 
while for k > k* the flow has a line of stagnation points. On the other hand, if h 
and k are fixed, then there exists a unique 7* > for which equality holds in (|5.19|) . 
which corresponds to a laminar flow stagnant at the bottom, and for 7 < 7* the flow 
does not have stagnation points, while for 7 > 7* the flow has a line of stagnation 
points. Whenever (|5.19|) holds, the stagnation line is Y = Yq, where 



tanh(A;/i) /tanh 2 (/c/i) g tanh(kh) 

(5 ' 20) h - Yo = ^k— + y—u^- V^T- 

from which we see that the distance between the stagnation line and the free surface 
decreases with k and with 7 2 , and approaches if and only if k — > 00 and 7 2 — > 00. 

Suppose now that 7 < 0. This case is the mirror image of the previous one. The 
flow corresponding to A_ never contains stagnation points. The flow corresponding to 
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A+ contains stagnation points if and only if (|5.19|) holds, in which case the stagnation 
line is Y = Yo, where Yo satisfies (|5.20p . 

It is of interest to compare Theorem 15.21 with the local bifurcation results for waves 
of constant vorticity 7 that were obtained in [T2] and in Let us fix k > 0. 

Varying h > in Theorem l5.2l we obtain all periodic traveling wave solutions of small 
amplitude and minimal period 27r//c that lie in a neighborhood of a laminar flow: 
for each h we obtain in the space (h, m, (//, w)) £ (0, 00) xlxX two smooth curves 
through the points (h,m-(h),0) and (h,m + (h),0) with m T given by (|5.14p . The 
union of these curves gives a surface M C (0, 00) xlxl 



X 




Figure 4 . The case Y > : The dashed curves represent the laminar bifurcating flows . 

Each nontrivial bifurcation curve ( in bold ) lies in a hyperplane h = constant . 

In Figure 4 we depict the situation in the case 7 > 0: 

(i) In the (h, m)-plane the curve h *— > m + (h) starts at the origin and approaches 

7/1 2 

the parabola m = as h — > 00, while m + (h) > for h > 0. As pointed out in 

the lines preceding (|5.19p . the nontrivial flows that bifurcate from the laminar flows 
corresponding to m+(h) do not present stagnation points. 

(ii) In the (/i,m)-plane the curve h 1— > m—(h) also starts at the origin and ap- 

7/1 2 

proaches the parabola m = — — as h — )■ 00. Denoting by h* > the unique positive 
solution of the equation 

tanh(kh) 7 2 h 

kh ~ 4g + 27 2 h 

we see from (|5.14p that m_(/i) < for h £ (Q,h*) and m_(/t) > for h > h*. 
Moreover, since 

7 2 K 7 2 h* 

±g + 2 7 2 K < g + 7 2 K 

1 1 r- r tanh(kh) . . . . . h 

and the function h 1— >• — is strictly decreasing on (0, 00) while h 1— > ^— - 

fC/i g + 7^ n 

is strictly increasing, we deduce that the nontrivial flows bifurcating from the laminar 
flows at m_(/i) contain stagnation points if and only if h > /i*, where h* G (0, /i*) is 
the unique solution to the equation obtained by requiring equality in (|5.19p . 

The approach in |44] is very different from ours, in that the author performs a 
/ Y \ 

flattening (X, Y) 1— > I X,— — - of the free surface Y = n(X). (This method is 
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therefore restricted to waves with non-overhanging profiles.) The bifurcation- inducing 
laminar flows identified there are exactly the same as ours. However, from the family 
of laminar flows of any depth ho, the local bifurcation curve constructed in [33] issuing 
from {hQ,m±(ho)) is a curve on A4 for which 



(5.21) Q = 2gh + 



m 7/i( 



2 



for the corresponding solutions of but no information on the conformal mean 

depth of the fluid domains is obtained (so that it need not be ho). There is a remark 
in [33] that a modification of the method there could have led to the existence of a 
new curve on M, for which, instead of (|5.2ip . the condition 

k f v / k 

(5.22) h o = ^ dx 

27T J-7T/k 

will hold for the corresponding solutions of . 

The approach in [T2], while more general in that it covers also non-constant vor- 
ticities, requires that there are no stagnation points in the flow. More precisely, it 
relies on a hodograph transform, for the existence of which it is essential that the 
horizontal velocity never vanishes in the flow. In [12] the relative mass flux m is fixed 
and h varies, so each nontrivial bifurcation curve on M is contained in a hyperplane 
m = constant. Since [12] only studies solutions for which i^y < in f2, it only obtains 
A_ in (|5.14p as a bifurcation point, and only in those cases for which the associated 
laminar flow (|5,15p does not have stagnation points. (However, the same method 
is easily seen to be applicable also for solutions for which ipy > in 0, leading to 
A + in (|5.14p as a bifurcation point, though again only in those cases for which the 
associated laminar flow (|5.15p does not have stagnation points.) In |12|. Section 3] it 
was concluded that, for k = 1, local bifurcation from trivial flows always occurs for 
negative constant vorticity, while in the case of positive constant vorticity there is a 
necessary and sufficient condition on the value of the mass flux m < 0, namely 



(5.23) tanhL/— ^) > 



7 / g + j^-2mj' 

Let us check these results against ours. We are thus looking at the laminar bifurcation- 
inducing flows for k = 1, and with tpy < everywhere in the closure of the fluid 
domain. Hence, by (|5.13p . (|5.14p and (|5.18p . and dropping the subscripts, 



(5.24) A = - - Jt^m + g tanh(/l) , 



(5.25) m = ?f - 7/lta 2 nhW - h ^ 72tan 4 h2(/t) + g tanh^), 

(5.26) \ + -yh<0. 

We are interested in the range of the mapping h \— > m(h) given by (|5.25p . restricted 
to the set of h £ (0, oo) for which (|5.26p holds. We will prove that the range is 
(—oo,0) if 7 < and coincides with the set for which f|5.23j) holds, thus confirming 
the conclusions of [12J. The case 7 < is easy, so we concentrate on the case 7 > 0. 
Recall that (|5,24p was obtained from (|5.12p with n = 1 and k = 1, hence 

(5.27) X 2 = (g - 7A) tanh(/i), 
as the solution with 

(5.28) A < 0, 
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while (|5.25|) was obtained from (|5.24|) by, see (|5.5a|) . 
(5.29) m = Ah+-±— . 

Hence m is in the required range if and only if (|5,27p and ()5.29p are solvable with 
respect to A G R and h > so that (|5.26p and (|5.28p are satisfied. We solve first 
()5.29p with respect to h for fixed A satisfying (|5.28p . The condition 



(5.30) A < -^-27™ 

is necessary and sufficient for solvability, and the requirement ()5.26p leads to the 
unique solution 



(5.31) /i = _A_v^+2^ ; 

7 7 

and to strict inequality in (|5.3U|) . Plugging (|5.31|) into (|5.27p we obtain 

A yj\ 2 + 2 7 m\ A 2 



(5.32) tanh 



7 7 / 5-7^' 



to be solved with respect to A in the interval (—00, —^/—2jm). The necessary and 
sufficient condition for this solvability is easily seen to be exactly (|5.23p . as required. 

Appendix A. On conformal mappings for L-periodic strip-like domains 

We prove here the existence and uniqueness of the conformal mean depth for any L- 
periodic strip-like domain For any such domain Q, let Q R be the domain obtained 
from £1 by symmetrization with respect to the real axis B, and let 

n = nuBun R . 

Let H > be such that £1 is a subset of the horizontal strip 

S H = {(x,y)eR 2 : -H < y < H}. 

The application 

-2-2 1 

ezH z — 1 

6 2H 2 + 1 

maps Sh conformally onto the unit disc D = {£ G C : |£| < 1}. Let II = <$(f2). The 
boundary of II is a Jordan curve J C B U {±1}, symmetric about the real axis and 
passing through the points ±1. Clearly $ is a homeomorphism between the closures 
of and of IT. Moreover, if D + denotes the part of B situated in the upper half-plane 
and n_|_ = II n B + , then II + = <I>(fi). Caratheodory's theorem |29t Theorem 2.6 and 
Corollary 2.6, p. 24] ensures the existence and uniqueness of a conformal map $1 
from B onto II, which has an extension as a homeomorphism between the closures 
of these domains and is such that $i(±l) = ±1 and <&i(0) = 0. The uniqueness 
assertion implies that ^i(^) = $i(z) for all z£D, This in turn implies that $1 maps 
the segment [—1,1] C M onto itself, and the imaginary part of $1 does not change 
sign in B + . Since the real part of <3?i is strictly increasing on the segment [—1, 1] C R, 
the Cauchy-Riemann equations on that segment imply that the imaginary part of $1 
is positive in B + . Hence $1 maps B + onto n + . Denote $0 = < &~ 1 - Since 

e f ( 2 +0 _ 1 

maps IZ2 conformally onto B , we see that $o°^ > i°^ > 2 = Uq+Wq with Uq, Vq : IZ2 — > M, 
is a conformal mapping of IZ2 onto O, with an extension as a homeomorphism between 
the closures of these domains, such that TZi is mapped onto and points symmetric 
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in 7^-2 with respect to the line y = — 1 are mapped into points in £1 that are symmetric 
with respect to the real axis. 

Let j3 + iS be the inverse of Uq + iVo, with f3, 5 : Q — )■ M. Since Q is an L-periodic 
strip-like domain, we observe that the function 

(X,Y) ^ S(X + L,Y) - 5(X,Y) 

is harmonic in f2, and has a continuous extension to the closure of which is zero 
on the boundary of 0, since by construction 5 = — 1 on B and 5 = on S. Being also 
bounded in £1 (since by construction |<5| < 1), the Phragmen-Lindelof principle [17] 
ensures that this function is identically zero, so that 

(A.l) 5{X + L,Y) = 5{X,Y), (X,Y)eU. 

The Cauchy-Riemann equations then imply that 

(A.2) 0{X + L,Y) - 0(X,Y) = K, (X,Y)eQ, 

for some constant K G R. Note that necessarily K > 0, since the fact that the real 
part of $i is a strictly increasing homeomorphism of the segment [—1, 1] CD implies 
that x i — y f3(x, 0) is a strictly increasing homeomorphism of M. From (jA.ip and ()A.2|) 
we infer that 

(u (x + K,y) = U (x,y) + L, 
\Vo(x + K,y) = V (x,y), 



(A.3) 



{x,y) G 111 



Replacing Uq + Wq from 1Z\ onto VLhy U + iV : IZh — > Q given by 

U(x,y) = U (x/h,y/h), V(x, y) = V (x/h, y/h), (x,y)eKh, 

where 

h = L/K, 

we see that h satisfies the definition of the conformal mean depth of 



p + i 8 



3> 



... 




U + iV 

o o 



Successive steps in the construction of the conformal map 



Suppose now that h also satisfies this definition, and let U + iV : — > £1 be the 
corresponding conformal mapping. Then 

(U + iVy 1 o (U + iV) 

is a conformal bijection from IZh onto 7?.^ and a homeomorphism between the closures 
of these domains. Let us write 

(U + iV)~ l o(U + iV) = t + iu, 
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where r, u) : IZh ~^ R are harmonic functions. It follows that 

(A.4) Ux + Ly)=r(x,y) + L, 

\u(x + L,y) =uj(x,y), 

Also, 

lu(x, 0) = 0, ui(x,—h) = —h, x £ R. 
The maximum principle shows that 

uj(x,y) = hy/h, (x,y)eK h , 

which implies that 

r(x,y) = (hx/fi) + C, (x,y)£U h , 

where C G R is a constant, a fact which, together with the first relation in (|A.4p 
implies that h = h. 

If, in addition, 5 is a curve of class C 1 '", then the Jordan curve J is Cj everywhere 
with exception of the two points ±1. By the Kellogg- Warschawski theorem [29} 
Theorems 3.5 and 3.6, p. 48], the conformal map $1 inherits this smoothness up to 
the boundary, namely it is of class C. in the closure of B with the exception of ±1 
and, moreover, its derivative is not zero at any point in D \ {±1}. Our construction 
ensures, taking into account the periodicity, that U, V 6 C l,a (1Zh) and 

U 2 X +V;V0 in W h . 

We now show that h coincides with the 'mean depth' of a L-periodic strip-like 
domain as defined in [TJ. There, the domain Q is considered by as being occupied by 
a fluid in steady irrotational motion, so that S and B are streamlines for a stream 
function ip which is L-periodic in X throughout Q and satisfies ()l.ld|) with 7 = 0, 
(jl.lep and (jl.lf|) . Since the flow is irrotational, there exists an associated velocity 
potential eft in $7, which is a harmonic function such that (j) + iip is holomorphic. 
The velocity field is then given by (ipy, — ifix)> which also equals by the 

Cauchy-Riemann equations. The mean horizontal velocity of the flow is then given 
by 

c = \ J <Px(T,Y)dT = ^ {4>{X + L,Y)-4>(X,Y)), 

a quantity which is easily seen not to depend on (X, Y) as long as the horizontal line 
segment joining (X,Y) and (X + L,Y) is contained in Q. The 'mean depth' is then 
defined as m/c, the ratio between the flux and the mean horizontal velocity. With 
(3, 8 as previously defined in Appendix A, note that the maximum principle implies 
that 

tp(X, Y) = m5{X, Y), {X, Y) € fl, 
and then the Cauchy-Riemann equations show that 

<f>(X + L,Y)- (f>{X, Y) = m(l3(X + L,Y) - f3{X, Y)), (X, Y) G 0. 
Therefore 

cL = mK 

and hence 

h = m/c. 

Consequently the conformal mean depth of £1 coincides with the 'mean depth' defined 
in p. 
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Appendix B. A commutator estimate 
Lemma B.l. If f G and g G C 2 ^* mt/i a G (0, 1), i/ien 
fC(g)-C(fg)eCi* forallSe(0,a). 

Proof. Let 

Q(t) = f{t)C(g)(t)-CVg)(t) foralliGR. 
Using the definition of C we obtain that, for all t G M, 

(B.l) 0(i) = i- jT cot {/(t) - /(,)} 5 ( S ) cfa, 

where, since / G C^, the integral is not a principal value integral, but a genuine 
Lebesgue integral. In proving that G C^f for all 8 G (0, a) we will make use of the 
following inequalities: 

1 i 

(B.2) \x - sin(x)| < -b| for all x G M, 

6 

(B.3) — ^— — < <t(|x|)— tt for all x with < Ixl < 7r, 
sin z (x) x z 



(B.4) 



/ \ 1 
cot(x) 

X 



< s(\x\)\x\ for all x with < Ixl < ir, 



where a,q : (0, it) — > [0, oo) are bounded functions on (0, a] for each a G (0, n). 

We show first that is differentiable on M and its derivative is given, for all i G 
by the formula obtained by formal differentiation under the integral sign: 

1_ r f'(t) S m(t-s) - {f(t) - f(s)} 
47r sin 2 



(B.5) 8(i) = — / _ _ 5(s)da. 



Note that the above integral finite. Indeed, the integral 

4vr J t _ T sm z (^) 

is finite, as it can be seen upon using (|B.2p . ()B.3p and the Mean Value Theorem, by 
which 

\f'm -s)- {f(t) - f(s)}\ = \t-s\ \f'(t) - 

< C\t-s\ 1+a , 

for some £ = s) with |£ — 1| < |s — t|. Hence the integral in (|B.5p is finite, being 
obtained from that in (jB.6[) by a rearrangement using the periodicity of the functions 
involved. In the above and throughout what follows, we denote by C a constant which 
depends only on the norm of / in and the norm of g in C 2 ^, and is allowed to 
vary from line to line. 

Since the operator defining commutes with translations, it suffices to show that 
is differentiable at t = 0, with derivative at t = given by (|B.5p , Note that, for 
any t close to 0, one can write 

where 

h{t) = \ f_l t { cot (nr) {/(*)-/(«)}- c °t (^) {/(o) -/(*)}} *00 (fa, 



-2|t| 

(-7r,-2|t|)U(2|t|,7r) 



/ f jcot ( ^ J {/(*) - /(,)} - COt L± ) {/(()) - /(.,)} \ n{ 
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Using (IB,4p . one can write I\(t) as the sum of a term which tends to as t — > and 



/_ 



2\t\ 
2\t\ 



[ f(t)-m f(o)-m 

\ t-s - s 



g(s)ds 



Since, for fixed t, the integrand in the above is continuous as a function of s, the 
Mean Value Theorem for integrals shows that the above quantity equals 

s 1*1 f f(t)-HO /(Q)-/(0 1 m 

for some £ G (— 2\t\, 2\t\). Using again the Mean Value Theorem, the above quantity 
equals 

for some ^1,^2 G (— 2|i|,2|t|). Since /' is continuous and g is bounded, the above 
quantity tends to as t — > 0. We have thus proved that 

(B.7) h(t)^0 asi^O. 

We now prove that, as t — > 0, 

/'(0)sin(-*)-{/(0) -/(*)} 



(B 



h(t) 



2sin 2 (f) 



g(s) ds. 



For this, we use the Dominated Convergence Theorem. Observe that l2(t) can be 
written as an integral over (— vr,7r), with the integrand given by 



X(— zr,-2|t|)u(2|t|,7r) 



/(*) - /(O) 



cot(M) + {/(0) -/(«)} 



cot(^) - cot(- f ) 



g(s) 



t K z ' " v ' ' v /J f 

where stands for the characteristic function of the set A. Firstly, note that, as 
t — > 0, the above integrand converges pointwise to the integrand on the right-hand 
side of (|B,8p . Secondly, note that, for each s 6 (— 7r, 0) U (0,7r) and for each i with 
< \t\ < \s\/2, the Mean Value Theorem yields that 



(B.9) 



<■"< I — — ) {/(*)-/(«)} 



cot 



{/(0) -/(*)}} 



/'(Osin^ 



{/(0 - /(*)} 



2sin 2 f£^ 



for some £ between and t. Using again (|B.2j) and ()B.3p . we obtain 



/'(0sin(£ -,)-{/(£) _/(,)} 



2 sin 2 



2 / t-s 



< C 



c 



/'(0(e --)-{/«) -/(«)> 



/'(0 - /'(&) 



+ C\£-s\ 



< cie-si"- 1 

< ci^r -1 , 



+ C|C-a| 



where £ is between and t, £0 is between £ and s, and we took into account that 
< \t\ < \s\/2 and / G C 2 ^\ Since s i-> |s|° _1 is integrable and 5 is bounded, we 
have thus checked the conditions of the Dominated Convergence Theorem. Therefore, 
(|B.8p holds and, since ()B.7p holds, it follows that O is differentiable at t = 0, with 
derivative given by ()B.5p evaluated at i = 0. As noted earlier, this implies that is 
differentiable on R, with derivative given by (|B.5p . 
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We now prove that 0' G C^f for each 5 G (0, a). Since 0' commutes with transla- 
tion, it suffices to show that 

(B.10) \Q'(t) - 9'(0)| < Cltl 5 for all t close to 0. 

The change of variables s t— > (t — s) in (|B.5p leads to the following formula 

mm fi'm 1 f gM^W - {/(*)-/(*-*)} w 
(B.ll) 9 (*) = - g(t - s) ds, 

for all t G R. It follows that, for each i close to 0, one can write 
(B.12) 4vr (e'(t) - 9'(0)) = Ji(t) + J 2 (t) + J 3 (t), 

where 

Ji{t)= . 2 / s \ g{t-s)ds 

J-2\t\ sm z (I) 

rm / 2| " -H ) sin(s) - {/(0) - f(-s)} 

M t ) = - . 2 fs\ g{-s)ds 

J-2\t\ sin^ (§) 

7'(*)sin(s)-{/(t)- /(t- S )} 
sin 2 (f ) 

/'(0)sin( S )-{/(0) -/(-,)} 



■w) = / \ ; 2~7T\ 9(t-s) 

./(-Tr,-2|l|)u(2|t|,7r) L sm {2) 



ds. 



J s {f(t)-f(t + T)}dT 



sin 2 (I) 

Note first that 
(B.13) 

r 2 ^ \sm(s)-s\ r 2 ^ 

\Ji(t)\< / \f'(t)\\g(t-s)\ l . V 1 + / \ 9 (t-s)\ — - a d S 

J-2|t| sm z (§) y_ 2 | t | sm z (| J 

<C|t| 2 + C7|t| Q 

where we have used (jB.3[) and (|B.2[) . the boundedness of /' and g, and the fact that 
/' G C^". Then, exactly the same argument gives that 

(B.14) I J 2 (t)\ < C\t\ a . 

Also, one can write 

(B.15) J 3 (t)=K 1 (t) + K 2 (t), 

where 

K(A [ {f(t) - rm sin( 8 ) + {/(o) - /(-,)} - {/(<) - /(t - ,)} 

K 1 {t)= r^y— g(i 

J(-tt ,-2|t|)U(2|t|,7r) sm {2) 

*M = / /-(0)sM») { /(0) -/(-.)} _ _ rfs 

J(-w -2|t|)U(2|t|,7r) sm z (2) 

Using ()B.3p . the boundedness of 5 and the fact that /' G C^, it follows that 
(B.16) 

I s (f'(r)-f'(t + T))dr\ 



I Isl Isl 2 



■/(-7r-2|t|)u(2|t|,7r) 

< C|*r I rids 

^(-7r,-2|t|)U(2|t|,7r) l s l 

<C\t\ a \ log |«|| 



(is 
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Also, using (|B.3p . the boundedness of /' and the the fact that g G C 2 ^, it follows 
that 

r r If (0)1 I- s f'( T ) dT ) 

(B.17) \K 2 (t)\<C | 5(t _ s )_ 5 (_ s) |/L^ + \ ds 

J(-tt ,-2|t|)U(2|t|,7r) L l S l l S l J 

< C\t\ a f ^- ds 

J{— Tr,-2|i[)U(2|t|,7r) \ s \ 

< C|t| a |log|t|| 

< C\t\ s . 

Combining the estimates (|B.12p - (|B.17p yields (|B.10|) . As noted earlier, this implies 
that 0' G • This completes the proof. 

□ 
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